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Preface 


The purpose of this memorandum is to discuss the evolution of wave energy, 
enstrophy and action for atmospheric Rossby waves in a variable mean flow. The 
presentation is theoretical, but does not represent original research; rather, 
it is pedagogic in nature. The work of a number of people has been drawn 
together into a unified account, with much of the algebra implicit in previous 
work made explicit here. The central results are the conservation laws given 
in Equations (4.24), (4.33), (4.34) and (5.10). 

It is a pleasure to thank David Andrews, who introduced the author to the 
concept of wave action and stimulated his interest in many of the concepts 
treated in this memorandum. Useful discussions with Brian Hoskins and 
Eugenia Kalnay are also acknowledged. Thanks also go to Ms. Mary Ann Wells, 
who did such an excellent job of typing the manuscript. 
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1. Introduction 


Conservation laws play an important role in geophysical fluid dynamics, 
for they yield strong insights into how flows are dynamically constrained. 

Such laws may apply to the fluid system as a whole, or to wave-like excitations 
of the fluid. The purpose of this article is to present a comprehensive 
development of the conservation laws of energy, ens trophy and wave action for 
atmospheric Rossby waves. This work is not new or original, but synthesizes 
and expands the published derviations of many others (Young and Rhines, 1980; 
Lighthill, .1978; Edmon et al ., 1980) into what is hopefully a clear and complete 
picture of how these conservation laws arise for Rossby waves, and what their 
limitations are. The powerful variational methods developed by Bretherton and 
Garrett (1968), Hayes (1970) and Andrews and McIntyre (1978) have not been used, 
as they are excessively abstract for the current purpose. 

To motivate the study of wave action and wave enstrophy, let us consider 
the more familiar notion of wave energy. Keeping in mind that we are dealing 
with linear systems, it is usually the case that wave energy in not conserved. 

For unstable waves, this is to be expected, since the terms representing the 
conversion of energy (kinetic or available potential) from the basic state to 
the wave will appear as a source term for the wave energy. However, conservation 
of wave energy also fails to hold for linear Rossby waves propagating in a 
slowly varing mean flow. There the local amplitude and energy density of the 
wave increase as the wave propagates towards larger values of the mean flow u, 
even though the wave does not draw energy from the mean flow (i.e. the wave is 
neutral). Since Rossby waves play a role in a variety of adjustment problems, 
including the mediation of tropical-extratropical interactions (see, for example, 
Hoskins and Karoly, 1981), it is of interest to find a wave quantity which _is 
conserved in a variable mean flow. The two central results of this article 
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are that: (1) In the absence of longitudinal (or x) variations in the mean 
flow, both wave action and wave enstrophy are conserved locally; and (2) For 
arbitrary slow variations in the mean flow, only wave enstrophy is conserved 
locally. The directional dependence of result (1) above is due to the fact 
that only meridional (or y) gradients of the planetary vorticity f are present; 
the anisotropic nature of the wave appears already in the fundamental dispersion 
relation. 

The basic strategy is to develop the elementary theory of a Rossby wave in 
a constant basic state u, and then to assume that even when u is variable, the 
wave dispersion and wave energy density are related to the local value of u 
and the local wave amplitude as in the elementary theory. The kinematic theory 
of Lighthill (1978) is then used to describe the rates of change of frequency 
and three dimensional wavenumber. This approach does not explicity make use 
of WKB theory (as does the derivation of Young and Rhines, 1980), although the 
spirit of the calculation is much the same. We also do not follow the approach 
of Hoskins e t al .(1983), who derive the local conservation law for wave potential 
enstrophy (our eq. (4.34)) because this approach does not make explicit how the 
frequency and wavenumber vector change in response to a slowly changing mean wind. 

It is perhaps important to point out that the local approach used in this 
article is not valid for unstable waves, since growth rates are generally 

determined by' the- conf iguration- of— the- mean - flow- u over- the whole domain- and -so 

depend on u in a non-local fashion. Thus, in order to study the form of the 
source terms that would arise (in the unstable case) in the conservation laws 
for wave action and wave enstrophy, one would have to either assume a vanishingly 
small growth rate, or seek an approach different than the one presented here. 

Section 2 presents the governing equations and the elementary theory of a 
Rossby wave in a constant "mean flow YiT TheT kinematica! - derivat icm“ for” the 
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rates of change of wavenumber and frequency for a slowly varying wave train is 


presented in Section 3. Using the results of Section 3, the elementary theory 
of Section 2 is generalized in Section 4 to the case of a slowly varying mean 
flow, leading to the derivation of the conservation laws for wave enstrophy and 
wave action. Wave action is related to wave activity and the Eliassen-Palm 
flux vector in Section 5. 
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2 • The Elementary Theory of Rossby Waves 
2a. The governing equations. 

The governing equations are those of quasigeostrophic theory on a beta 
plane. They are most easily derived from scaling considerations when altitude 
(z*) rather than pressure (p) is used as the vertical coordinate (e.g. Charney 
and Stern, 1962). Furthermore, Rossby waves are more easily treated using z* 
as a vertical coordinate. On the other hand, the advantages of using p as a 
vertical coordinate are well known. In this article we follow Palmer (1982), 
and retain the advantages of both systems by using log (pressure) coordinates: 


z = -H £n (p/p Q ), (2.1) 

where p Q = 1000 mb is a reference pressure and the scale height H is a constant. 
The quasigeostrophic equations on a beta plane are derived in z coordinates in 
the Appendix. They take the form of the conservation of potential vorticity 
following the horizontal motion of the flow: 

(JL - 14 JL + 14 A.) q = 0 , (2.2) 

9t 9y 9x 9x 9y 

where the potential vorticity q is given by: 

q = f + 3y + V% + f ° H e Z/|,H JL (e _z / H ii) (2.3) 

° 2 RS 9z 9z 


Here is the geostrophic streamfunction, defined so that u = - v = 

<5y dx 

2 

with u being the eastward and v the northward velocity, and is the horizontal 

9 9 

Laplacian operator f— — + ~ — 1. Also, f = 2Q sin*, where Q is the angular 

9x 2 9y 2 

rotation rate of the earth and <b • the central latitude of the beta plane, 6 = — 

° dy 

where y is defined by dy = a d(j), a being the earth* s radius, f = 2ftein<J>, R is the 

<v 

gas constant, and S = (£— ) K ^s. k = R/C p , C p being the specific heat at 


V 
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constant pressure, and 0 S = 0 s (z) is the reference state potential temperature 
(Charney 1947: Charney and Stern, 1962). In deriving (2.3) it has been assumed 
that S is constant in the vertical. We shall make that assumption throughout 
this article. 

f 2 H 

It is useful to relate the parameter ° appearing in (2.3) to the Brunt- 

RS 

Vaisala frequency N_, given by N ^ = i£— dJs where z* is altitude, as before. 

s s e s dz* 

Since ills, fL®s and d 9 S depend only on the reference state, 
dz* dp dz 


Ik - Ik - 

dz* d P dz* 


P s g = 


_d 0 S pg 

dp” R?7 


where d 8 S (written with p as an independent coordinate) should be written as 
dp 

dQg if z* is the independent coordinate. Similarly p apearing in the above 

dps 

formula should be understood as p s in the z* system. We have used the equation 

2 

of state of an ideal gas, p = pRT. From this, N„^ = .IEJL — 4j*s. Introducing H_ 

S RT S 0 S dp s 

the scale height of an isothermal atmosphere with temperature T s 


"Pg d0 g gH d0 g gH 

®s^s dp ®s^s dz H s T £ 


(2.4) 


so that S - -5-i N 2 
gH S 


and 


RH T N 

RS = s s s 


f o H 8 f o H ! 



(2.5) 


where equation (2.5) defines S. This quantity will appear in many of the 
subsequent equations in this article, and is closly analagous to the parameter 
S introduced by Pedlosky (1979, p. 334), who developed the equations in z* 
coordinates • 

2b. Linearization: The Prototype Rossby wave 

Equation (2.2) is linearized about a basic state (denoted by an overbar): 
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q = q + q ' , and 

^ . 

Here q and ^ can be functions of (x, y, z) but not time t, whereas the 
perturbations q^ and ty' are functions of (x, y, z, t). Note that 

? = f Q + gy + v| t + 1 e z/H JL (e _z/H il) , (2.6a) 

q " = V? r + A e z/H A_ (e -z/H ill) . (2.6b) 

3z 3z 

The most general linear form of (2.2) is 


j q ^ + u ^q + u ^ + v + v " J9. = 0. (2.7) 

3t 3x 3x 3y 3y 


We will return to (2.7) in Section 4. The prototype Rossby wave is obtained in 
this section by setting 

- = -uy, (2.8) 

where u is assumed constant. Thus v = 0, and from (2.6a) q = f Q + 0y, so that 

_5l = 0 and ** 6 . Equation (2.7) becomes 
3x 3y 


hi + u hi + = 0 . (2.9) 

3t 3x 


We “try a solution of the form 


ty' - Ae z /(2H) cos (k x + ly + mz - wt) = Ae z /( 2H) c os( a) . (2.10) 


Here A is the (constant) wave amplitude, (k, £, m) the (constant) wave vector, 

and a) the (constant) frequency. 

From (2.6b) and (2.10) it is easy to show that 
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= -[k 2 + £ 2 + ( m 2 + _1 _) ] Ae z /(2H) cos ( a ) 

^ 4H 2 

= -K 2 Ae z /( 2 ^) cos( a) , (2.11) 

which defines the three dimensional wavenumber 

K 2 = k 2 + £ 2 + 1 [ m 2 + J_] . (2.12) 

^ 4H2 

Using (2.10) and (2.11) to evaluate each term in (2.9), 

ill = - 0 ) K 2 Ae z /(2H) sin( a) , 

3t 

ill * +k K 2 Ae z /( 2H > sin(ci) , 

3x 

and v' = ^ = -k Ae z /( 2 ^) sin( a) • 

3x 

Equation (2.11) then gives the dispersion relation of a Rossby wave propagating 
in a constant mean flow: 

a) = u k - JL (2.13) 

K 2 

2c. Energy and Potential Enstrophy Density 

The perturbation energy is defined as the sum of perturbation kinetic plus 
available potential energy. In ordinary pressure coordinates (Lorenz, 1960), 

e tot = e k + e a > 

with e k . ;//d*dy &. i [(|C) 2 + (|t) 2 ] 

g 2 9x 3y 

and Ea = // /dxdy ^ , 

S 2S 

_ ,D x K d0 

where S = -p ( J s = HS . 

Po d P 
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Converting from p to z coordinates via (2.1), 


e tot - % ///*>■ »y ( i l(|f f * (| f ) 2 ] + \ ^P- 2 } • 

gH 2 dx dy 2 s 

f H 

Using the quasigeostrophic relationship T' = 0 ° v 1 ' (see equation (A. 10) in 

R 3z 

the appendix), and introducing T = JlS. and p. i — = _J-2. , we have 

0 R ° gH RT 0 

E I0T - // /dxdy d, p 0 .-/« ( I '[(iff + (if) 2 ] + i (it) 2 ] 

2 3x 3y ^ 3z (2.14) 

The temperature T 0 is that which renders gH/R a scale height for an isothermal 
atmosphere, while p 0 is an "equivalent" surface density. 

From (2.14) we define an energy density (energy per unit mass) e as 

e = I e-z/H [{Hlf + [Hlf + I (ill) 2 ] . . (2.15) 

2 3x v 3y /g\ 3z J J 

Note that the energy density includes the exponential decrease of density with 
height. From (2.10), 

e = i. A 2 {(k 2 +£ 2 ) sin 2 (a) + JL [ 1 cos^(a) + sin2(a) 

2 ^ 4H2 

- iH cos(a)sin(a) ]} 

H 

Averaging over one wave period 


2tt 

— E ='_L~ / da e = I~A 2 (k 2 +Jl2 + -l [ _JL..+ m 2 - -])- = I A 2 K 2 

2n o ^ 't 4H 2 * 


(2.16). 


The total perturbation potential enstrophy is 


P TOT ~ ///dxd y ^2- (i q ) • 


(2.17) 


As with energy, this can be written as P^qT 


J/fdxdydz p Q e ~ z ^ H (l. q"^) , leading 
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to the definition of potential enstrophy density p as 


p = e- z / H I q' 2 . (2.18) 

. 2 

Using (2.11), and again averaging over one wave period, 

2n . 2 

P = i / d a p = J_ [ Acos ( a) ] , 

2tt o 2 




3 . The Kinematic Theory of Wavetrains 


Here we derive the fundamental properties of wavetrains, or almost purely 
sinusoidal waveforms, in the presence of constant or slowly changing mean winds. 
The material presented in this section is condensed from Lighthill (1978), and 
is necessary to the derivation of the conservation laws in the following section. 
The kinematic theory demonstrates the key role played by the group velocity, and 
is the basis of ray tracing, a technique which has recently found widespread 
application (Karoly and Hoskins, 1982). 

3a. Constant Mean Wind 

We let the mean wind have components in both horizontal directions, so 
that both u and v are non-zero, although they are both constant. In this case 
(2.7) reduces to 

dq I + u dq 1 + v ill + gv^ = o , 

3t 3x 3y 

for which (2.10) is a valid solution, with 

oj = (uk + v£) - _$S 
K2 

Here we do not need the specific form of the dispersion relation, we only need 
the fact that there is one, i.e. 

a) = ft (k, £, ra) . (3.1) 

In the presence of the constant mean wind (u, v), a more general solution 
of the linearized equations can be constructed as a superposition of solutions 
of the form of (2.10). The approach adopted here, however, is somewhat different. 
We consider a solution which locally is given by (2.10), but in which the amplitude 
(A), frequency ( u>) and vector wavenumber (k, £, m) change slowly in space and time. 
Physically, this corresponds to the assumption that one has waited long enough 
after an initial disturbance for the different components of the spectrum (waves 
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of different wavenumbers) to be separated in space from each other, i.e. for 
the disturbance to have dispersed. The amplitude, frequency and wavenumber 
components change by a small fraction of themselves over the local wavelength 
or local period. 

We start with the generalized form of (2.10) 

V - A (x, y, z, t) e z /(2H) cos ( a ( x , y, z, t)), (3.2) 

and expand the total phase a in deviations (6x, 6y, 6z, 6t) about some point 

( x o> yo> z o> t o) : 

a(x, y, z, t) = a(x , y , z , t ) + 12 6t + — 6x + 6y + 6z 

° ° 3t 9x 3y 3z 


+ 0 ((6t) 2 , ( 6x) 2 , ) + (3.3) 

We can, without loss of generality, set ot(x 0 , y Q , z Q , t Q ) = 0, and, 

neglecting second order terms: 

ot(x, y, z, t) = <5t + .^2 6x + 6y + ..^2 ^z . (3.4) 

3t 3x 3y 3z 

Since a locally takes the form 

a = kx + £y + mz - cat , 

the local frequency oj and wavenumber vector (k, ra) are defined as: 

w = - i« (3.5a) 

3t 

(k, l, m) = fi«, -1®, iS) • (3.5b) 

3x 3y 3z 

We now come to one of the most important ingredients in the theory, the 
step at which the physics of the wave enters the formalism. I t is assumed that 
the dispersion relation (3.1) holds for all space and time , with the frequency 
and wavenumber being given by their local values. From (3.1) and (3.5a), (3.5b), 
this gives: 
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( 3 . 6 ) 


-3a = 0 ( 3a , 3a ? 9a) 

9t 9x 9y 9z 

Differentiating (3.6) with respect to x, y and z, separately, gives: 


9 2 ot = 



+ 9fi J 

9 2 a ) 

^ 9Q | 

[i*.) , 

9x9t 

9k 

9x2 

9£ 

9x 9y 

9m 

9x 9z 

.2 

9 a . - 


(9 2 a 

) + 15 


+ 15 ( 

:ii«_) , 

9y 9t 

9k 

9x9y 

9£ 

9y 2 

9m 

9y 9z 

9 2 a = 

15 

(9 2 a ' 

) + 15 

( 9 2 a 

) + 15 

(ii«) . 

9z 9t 

9k 

9x9z 

9£ 

9y 9z 

9m 

9z2 



-> 






If we define the group velocity C as: 



C , C 
y ’ z 


) . (is, 

9k 


dn 

9£* 


9fl ) 

9m 


(3.7a) 

(3.7b) 

(3.7c) 


(3.8) 


and use (3.5a) and (3.5b) to derive relationships of the type 


9 a — _ 3 (jj — 9k 
9x9t 9x 9t 


9 « = . e tc. 

9x2 9x 


( 3 . 9 ) 


equations (3.7a) - (3.7c) become 


-> -> 


— + C • Vo k 
9t 


= 0 


1A 

9t 


9m 

9t 


-> -> 


+ C 


V 3 £ = 0 


-> -> 

+ C • Vo m = 0 . 


-> 


(3.10a) 


(3.10b) 


(3.10c) 


Here V 3 is the three-dimensional gradient operator. Thus, taking the total 
time derivative going with the group velocity all the wavenumber components 
remain constant. Invoking the definition 


D _ 9 


-> -> 
+ C • Vo 


-Dt 9t J 


(3.11) 
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here and throughout this article, (3.10a) - (3.10c) are just 


Dk = o 

(3.12a) 

Dt 


D* = 0 

(3.12b) 

Dt 


Dm = 0 . 

(3.12c) 

Dt 



Further, w is also constant following the group velocity, for 


D co = 9Q Dk + 30D&+ 90 Dm — q 
D t 9k Dt 9£ Dt 9m Dt 


(3.13) 


3b. Spatially varying mean wind. 

We are now in a position to introduce slow spatial variations in the mean 
wind, slow in the sense that, as before, there is a small fractional change 
over a wavelength. We assume the mean wind to be constant in time. The 
dispersion relation is still assumed to hold for all space and time, but since 
the dispersion relation generally involves the mean wind, which is no longer 
constant, (3.1) must be revised. The frequency w now has a (slow) dependence 
on the coordinates beyond that given by the (slow) coordinate dependence of 
(k, £, m). This additional dependence of u) on (x, y, z) is formally introduced 
by letting u) depend on the parmaters X^ , i=l, ...... N. (3.1) then becomes 

a) - SI (k, £, m; X ± ) . (3.14) 

For most applications, the X^ can be taken to be the components of the 
mean wind: \\ = u, X2 = v. Equation (3.6) becomes 

- is = a (is, is, — ; A, ) , (3.15) 

at 9x 9y dz a 

and (3.7a) - (3.7c) acquire the additional terms 
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-3a , 

. . .+ y 

30 

3 A-; 

(3.16a) 

9x 9t 

L 

i 

3Xi 

9x 

- 3 2 « = .. 

. . .+ y 

da 

3Xj 

(3.16b) 

3y 3t 

L • 

i 

3Xi 

3y 

_ 3 2 a = 

. . .+ y 

30 

3X.| 

(3.16c) 


3z3t i 3A^ 3z 


while (3.12a) - (3.12c) become 


(3.17a) 
(3.17b) 
(3.17c) 

Finally (3.13) becomes 

Dm = j)0 Dk + _30 D_£ + J)0 Dm + Y 30 D Aj # (3.18) 

Dt "3k Dt "3£ 5F 3m Dt i 3^ Dt 

From (3.17a) - (3.17c) and (3.18), we prove a remarkable result: 

Da) = _ 30 [V 30 3 X-; ] _ 30 [ £ J0_ _3 A* ] _ 30 [ £ _30_ 3X t ] 

Dt 3fc i 3Xj 3x 3£ j[ 3X^ 3y 3m ^ 3X^ 3z 


Dk = 

- i 

dtt 

3Xi 

Dt 

i 

3X ± 

3x 

D£ = 

- y 

30 

3Xj 

Dt 

L 

i 

3X ± 

3y 

Dm - 

- y 

30 

3 Aj 

Dt 

Lt 

i 

3Xi' 

3z 


+ y*L 

i 3Xi 



iFi + c 

9x 


y 



+ c„ 



9 


where we have used (3 #11). Noting that = 0 (mean winds constant in time), 

dt 

and "using (3*8), we easily find that 

Ow = 0 . (3.19) 

Dt 
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The most important source of dependence of Q on the arises from the 
advective effect of the mean flow. This dependence can be expressed by 
introducing the transformation to the frame of reference moving with the local 
mean flow , i.e. the frame of reference in which the local mean flow vanishes. 

The coordinates in the transformed frame are denoted by primes. We have 

x' = x - ut 

y" = y - vt 

(3.20) 

z' - z 

t' = t 

We obtain the transformed wavenumber (k' > V , roO and frequency (cO by noting 
that the local phase a is invariant, so that to first order in the derivatives of a, 

(- + k^x^ + Z' y' + xb'z') * (-at + kx + Zy + mz), (3.21) 

for all (x, y, z, t). From (3.20), this becomes 

[- (a + k^u + Z'v) t + k ^x + Z'y + m^z] 

= [- a£ + kx + Zy + mz] ,* (3.22) 

whence k = k ^ , Z - Z' y m = m x , and 

a) = o)^ + ku + Zv (3.23) 

a)" is called the intrinsic frequency. 

We may not have transformed away all the dependence on the variable mean 
wind, however, for it is possible that to" still depends on u, v. An example of 
this arises in the study of Rossby waves in shallow water theory (Pedlosky, 1979). 

In that problem depends on u via a term arising from the perturbation advection 
of mean state potential vorticity (Pedlosky, pp 91, 109). Physically, the mean 
current u is balanced by a slope in the y direction of the free surface, a slope 
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which is the same in any frame of reference moving in the x direction. In the 
case treated in the previous and subsequent sections, namely Rossby waves in a 
stratified fluid, is independent of (u, v), and depends on (x, y, z, t) only 
through (k, £, m). 

Analogously to (3.14), we can write 

(x)" = to' (k, £, ra; X t ) , (3.24) 

where again refers to u and v for i=l and 2. 

From (3.23) and (3.24), _ __ 

a) = ku + £v + 9,' (k, £, m; X^) (3.25) 


The right hand sides of (3.16a) - (3.16c) can now be written as: 


i 3 X-£ 3x 3x 3x 


y 3Xj 

i 9y 


y iiL iAi 

i 3Xi 3z 


3y 


k 9u + £ 3v + l 
3z 3z i 


3cr 

3 A{ 

(3.26a) 

axi 

3x 


3JT 

IXi 

(3.26b) 

3Xi 

3y 


3JT 

3X-t 

(3.26c) 

3Xt 

3z 



Here (k, £, m) have not been differentiated, since it is only the coordinate 
dependence beyond the dependence through (k, £, m) that is being considered. Finally, 
combining (3.26a) - (3.26c) with (3.17a) - (3.17c) and specializing to the 
case of Rossby waves in a stratified fluid f 9ft = 0 1, we have 

(3.27a) 
(3.27b) 
(3.27c) 


Dk 

= _k 3ii - £ 3v 

Dt 

3x 3x 

D % 

1 

1 3 

CO 

1 

II 

Dt 

3y 3y 

Dm 

= -k - l . 

Dt 

3z 3z 
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4 . Conservation Laws for Wave Energy, Wave Action and Wave Enstrophy in a 


Slowly Varying Mean Flow 
4a. Conservation law for wave energy 

We start with a full linearized equation (2.7), in which u and v may vary 
slowly in space: 


jSL + u dg + v 3q + u ^ +v * j|£L = 0 , (4.1) 

9t 9x 9y 9x 9y 


where q is given by (2.6a). With u =» - v = JL4, 

9y 9x 


8q = V 2 v + 1 e z / H JL (e~z/ H A.) (4.2a) 

3x 2 f 3z ' 3z 


and 

AL = B + V 2 u+1 e z / R JL (e -z / H Ail (4.2b) 

9y 2 ^ 9z 9z 


Since the mean flow is assumed to vary slowly in space, we neglect second 
derivatives of (u, v) in (4.1). 

Thus , 

Al - o, \ 

9x j 

Al * 0 , ) 

3y 


(4.3) 


and (4.1) becomes 

(A_ + u A_ + v A_1 q' + = 0 , (4.4) 

9t 9x 3y 


where q' is given by (2.6b). We try a solution of the form (2.10), but in 
which the amplitude A and vector wavenumber (k, m) are allowed to vary 
slowly in space and time . The variations of (k, £, m) obey the relations 
derived in Section 3, especially equations (3.9) and (3.27). The frequency to 
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is given by the local dispersion relationship, while the variation of the 
amplitude A is discussed is the following paragraphs* 

Our basic solution is 

= A( x , y , z , t ) cos(kx + iy + mz - cat) = A(x,y,z,t) e z /(2H) cos(a). (4.5) 

We substitute (4.5) into (4.4), using (2.6b). A number of intermediate results 
are helpful: 

= §A e z /(2H) cos(a) - kA e z /(2H) sin(a) (4.6a) 

8x 3x 

9 = 3__A e z/(2H) cos ( a) - 2k JA e z /(2H) sin(a) 

8x2 3x2 8x 


- A e z /(2H) s in( a) - k 2 A e z /(2H) cos ( a ) . (4.6b) 

8x 


A result analogous to (4.6b) holds for JLJLr, so that 

8y 2 

2 2 

V = V A e z /(2H) cos(a) - 2(k + £• — ) e z /(2H) sin(a) 

2 2 _ ^ 9x 9y 

- A C §j£ + JA) e z /(2 H) sin( a) - (k^ + £.2) ^ e z/(2H) cos(a) . (4.6c) 

9x 9y 




Also, 


1 e z/H j3_ f e -z/H 3^ ) = j_ e z/H [ e -z/H f3A cos ( a ) 

^ 3z 3z ^ 3z 3z 

+ A_ cos( a) - mAsin( a) ) 1 
2H 


= A e 2 /^ [- i e z / (2H) (3A cos ( a ) + A__ cos ( a) - A sin( a) ) 

^ 2H 3z 2H 


2 

+ Q-z/iZW) 1 9 A c os(a) - 2m JA sin(a) + J JA cos(a) - ®A sin(a) 

g z 2 3z 2H 3z 2H 

“AE A sin(a) - m^A cos(a) }] 

3z 

= J_ e z/ ( 2H) J cos(oi) • {AA - 1 A. — A + 1 9A _ m 2^ } 

^ 3 z 2 2H 3z 4jj2 2H 3z 


+ sin(a) • { lA - 2m JA - HA - A AE } 1 
2H 3z 2H 3z 


2 

= A [ cos(a) {1_A - A(m^ + A ) } + sin(a) { “2m j!A - A AE }] 

^ 4r2 3z < 4 - 6d ) 

Further, 

(A_ + u JL_ + v JL.1 cos(a) = ( a) - ku - £v) sin(a) \ 

3t 3x 3y J 

[ (4.6e) 

(A_ + u A_ + v JLA sin( a) = -( a> - ku - Zv) cos(a) / 

3t 3x 3y 


Since the amplitude is assumed to vary slowly, we will neglect the under- 


lined terms in (4,6b), (4,6c) and (4,6d), all of which involve second spatial 

2 

derivatives of A. For example, the neglect of j A in (4.6b) compared to 2k 

2 7 * 

is equivalent to the statement that AA A « JA, where A is the wavelength in 

3 X 2 3x 

the x direction, which is just the statement that the change of JA over one 

3x 

wavelength is small. 
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From (4.6a)-(4.6e), and (2.6b), 


q' = e z ! (2H) [cos(a) • {-A(k 2 + £2 + J_ ( m 2 + 1 ) ) } 

^ 4H2 


+ sin(a) { -A (ill + HL + 1 i®) - 2 (k^ + £ j>A + 1 ra JA) } ] . 

3x 3y ^ 3z 3x 3y ^ 8z (4.6f) 


Finally (4.4) becomes 

( (0 - ku - £v) • sin( a) • (-AK 2 ) + ( w - ku - £v) • cos( a) * [a fj?L + JL?: + }_ 1 

3x 3y IT 3z 

+ 2 fk_^ + £ i*A + _L m 9A 1 1 

3x 3y ^ 3z 


+ cos(a) [- K 2 dA - 2 A (k dk + £ d£ + 1_ m dm)] 
dt dt dt Q dt 


+ sin(a) [- dA ( 3k + _3£ + 1 3m j ^ d_ f 3k + _3£ + 3m j 

dt 3x 3y ^ 3z dt 3x 3y ^ 3z 


- 2 (dk_9A + d£^A + ldm _3A) _ 2 (k d \ 3A + £ d_ _3A 
dt 3x dt 3y ^ dt 3z dt 3x dt 3y 


+ JL m d_ _3A)] * + 3 j^A cos(a) - 6k A sin(a) = 0 . (4.7) 

^ dt 3z 3x 


In (4.7), K 2 is given by (2.12), 

K 2 = k 2 + £ 2 + i (m 2 + __L_) , 

^ 4H 2 


and d_ i s defined by 
dt 


l=i+ul. + vi . (4.8) 

dt 3t 3x 3y 


Both the operator d__ defined above and the operator J-L. defined in (3.11) will be 

dt Dt 

needed in the subsequent development. The term with the asterisk in (4.7) can 


be neglected, since it contains second derivatives of the slowly varying quantities 


- 20 - 



A, (k, £, m) or products of their first derivatives, (4,7) is really two equations, 
one for the coefficients of the sin(a) term, the other for the coefficients of 
the cos(a) terra. Setting the sum of all coefficients of sin( a) to zero, we get 

- (w - ku - £v) A K 2 - 8kA = 0 (4.9a) 


or 



(4.9b) 


which is just the generalization of (2.13) to the case of non-zero mean flow 
in the y direction. 

The equation obtained by setting the sum of the coefficients of cos(a) to 
zero in (4.7) is: 


( 0) - ku - £v ) [a (Ji + + 2. i®) + 2 (kiA + £ _9A + i m JA)] 

9x 9y / < 5 V 9z 9x 9y 'g' 9z 


-K 2 dA-A±_K 2 +8iA=o . 
dt dt 3x 

With the help of (4.9b), (4.10) becomes 


-gk ^ f 8k + d_Z + 1_ dm j + 2 
^2 3x 3y ^ dz 


(k + JL in 3 1 

3x 3y 'g' 3z 


_ k 2 dA . 

■ A L K 2 + 

8 JA = o 

dt 

dt 

3x 


(4.10) 


(4.1 la) 


or 

dA + 2Jk 3A + £ W + 1 m 3A) . £ 3A 

dt ^4 3x 3y 'g' 3 z k 2 3x 


= ZA iL_ - g kA [ 3k + dJL + J_ dm J 
](2 dt ^4 3x ^ 


(4.11b) 


-> 


To make further progress, we must calculate the group velocity C where, as 
in (3.8), 
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We have, from (4.9b) 


-> 


c = (c , c c ) = (™, iSJ, 1“). . 
x y x 3k 3£ 3m 


C = — = u - -i- + = u + C' , 

X 8k K 2 K* 


(4.12a) 


C v = — = v + = v + C ' 

y 3£ ^4 ^ 


(4.12b) 


^ _ 3w = + 2Pkm - p 

z ' 3m " ^4 ' 2 


(4.12c) 


-> 

Equations (4.12) define the intrinsic group velocity C but it is also given 
-> 

by C " = ( do) , Ail , AO , where us' = a) - ku - Iv is the intrinsic frequency as 
3k 3£ 3m 

defined in (3.23)- Using (4*12a)-(4.12c) in (4.11b), we get: 


+ c — + c , — + c, 9A 


3t 


3x 


3y 


3z 


DA = _ A_ d_ ( K 2) _ JkA [3k + ^ £ + 1_ 3ml . 

Dt £2 dt ^4 3x 3y / g v 3z 


(4.13) 


Now, from (2.16), 


DE=J_ AK 2DA + J_a 2D_ ( K 2) 
Dt 2 Dt 4 Dt 


(4.14) 


Multiplying (4.13) by i AK 2 , and adding _L a2 D_ (r 2) to both sides, we get: 

2 ' 4 'Dt - - - - ----- - 


DE = ^ A 2 D_ ( K 2) _ 1 a 2 d ( K 2) _ j_ gkA [J3k + dJL + 1_ jta] 

Dt 4 Dt 2 dt 2 j(2 3y 3z 

V J\ _ J 


(4.15) 


© 


© 


where it must be kept in mind that the operators iL_ and IL_ are not identical. 

dt Dt 
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Term 1 in (4.15) can be further developed: 


fl)= I A 2 — (K 2 ) - I A 2 — (K 2 ) + I A 2 (C' — + C' + C ' K 2 

4 Dt 2 Dt 2 x 3x y 9y z 3z 


= - 1 A 2 £_ ( K 2) + I A 2 [f Zi + iiii] .L K 2 + IM 1_ K 2 
4 Dt 2 

+ 2 ^ m 1_ K 2 1 , 

/ §K.4 9z 

so that 


(T)= -I A 2 [ k DJi + £ M + i m Dm ] + 1 k 2 [ f - fi + 2fr*_) 3_ K 2 
v 2 Dt Dt ^ Dt 2 ^2 ^4 £x 

+ 2gk£ i_ K 2 + 2 ^- m JL K 2 ] . (4.16) 

k 4 3y 'sk^ 9z 

Using (3.27a)-(3.27c) in (4,16), and putting the result back into (4.13) gives us: 


DE = 1 A 2 [k fk in + Z ill + Jtfk iH.+ £ ill + i m (k i“ + £ ill] 
Dt 2 3x 3x 3y 3y ^ 3z 3z 


+ i A 2 [fci + 2jd) 1_ K 2 + 2 ^ £ i_ K 2 + i 2 gkm K 2 ] 
2 ^2 j^4 8x ^4 dy ^ ^4 3z 


_ _1 Sc A 2 [ 3k + _3£ + .1 3m 1 (4.17) 

2 k2 3x 3y ^ 3z 

Equation (4.17) is more useful than (4.15) because there are no time derivatives 

on the right hand side. It turns out that the right hand side of (4.17) is 

-> 

related to the divergence of the group velocity C . Starting from (4.12a)-(4.12c) , 


”> “> 3C 3C 3C 

7_ . C = -^L + -JL + —L = 

J 3x 3y 3z 


3u + i_ (K 2 ) - 4 gk 2 3 (K 2 ) + 2_S _3_ (k 2 )+ - 4 gk £ 3 (k 2 ) 

9x ^4 dx ^6 dx j(4 3x dy %6 dy 


+ li JL (k£) - i_ (K 2 ) + 2JL 1- (km) . (4.18) 

k 4 3y 'sk 6 9z 9z 
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Multiplying the LHS of (4.18) by E and the RHS by its. equivalent , via (2.16), 


e Vo *c = — a 2 k 2 (in. + lx.) + I JL (k 2 ) - JM! [k JL (k 2 ) 

^4 3x 3y^ 4 ^2 3 x ^4 

o 

+ l JL (K 2 ) + 1 m JL (K 2 ) 1 + I [A. (k 2 ) + JL (k£) +11. (km) ] . (4.19) 

3y ^ 3z ^2 K 2 3x 3y '§'■ 3z j 

V © 

We now assume that the mean flow is non-divergent : 


JHi + _3v - q 

3x 3y 


(4.20) 


Then terro(^T)in (4.19) vanishes. Terra(2)can be simplified using identities relating 
the slow changes of wavenumber components which are derived from the fundamental 
definition of wavenumber given in (3.5b), namely 


d (k £) = k + £ ^k - k JL?: + £ ^ a ~ k ^ ^ ^ 

3y 3y 3y 3y 3y 3x 3y 3x 

and 

2 

j (km) = k + m ^k = k _3m + m 3 a = k ^m + m 
3z 3z 3z 3z 3z 3x 3z 3x 


(4.21a) 


(4.21b) 


With (4.21a) and (4.21b), terra(?)in (4.19) becomes 

2\ = 1 - BA 2 k [Jk ..+-34 + J_ 3m] +. JA 2 . [ k _3k + + I tn i® ] - - , - - 

^ 2 yA 9x 3y ^ 3z 2K^ 3x 3x ^ 3x 

= k [Jk + 11 + A 3m] + JA 2 ^ _3_ ( K 2) . (4.22) 

2^2 3x 3y 3z 4 j^ 2 3x 

Using (4.22), (4.19) finally becomes 
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-> -> 

E V 3 *C = 


BA 2 V [3k + U + 1_ 3m] + BA^ a_ 

2^2 3x 3y ^ 3z 2K^ 9x 

- gk A 2 [k JL (K 2 ) + £ JL (K 2 ) + i m i_ (K 2 ) ] , 

K 4 Sx 3y ^92 (4.23) 


which can be used to rewrite ( 4 • 17 ) as 

— = -E V, *C + -La 2 [k(k — + £ il) + £fk + £ il) 

Dt 2 3x 3x 3y By 

+ I m fk _§JL + £il]] . (4.24a) 

Q 3z 3z 

Remembering the definition of ]L_ (given in (3.11)), an equivalent form is 

Dt 


9E + 

-> -> 

V(CE) = 

= 1a 2 

[kfk i" + 

£ ix + £ 

0 + 

£ i£) 

at 


2 

ax 

3x 

ay 

9y 





+ _L m 

(k^ + 

£ iE)] 





'S' 

dz 

a z 


(4.24b) 

This is one of the important results in this article. It states that in the 
absence of mean flow gradients, the wave energy density satisfies a local 
conservation law giving the rate of change of energy in an arbitrary volume in 
balance with the flow of energy into or out of that volume, with the appropriate 
flow velocity being the group velocity . 

In the presence of mean flow variations, wave energy is not conserved . We 
must turn to wave action and wave enstrophy to find conservation laws which 
hold in the presence of mean flow variations. 

4b. Conservation laws for wave action and wave enstrophy 

It turns out to be convenient to define a generalized wave action by 


- 25 - 




(4.25) 


£= fOOj , 


where f(k) is an arbitrary function of the wavenumber component k, and the 

intrinsic frequency to' is given (by (3.23) and (4.9b)) by 

to' = , (4.26) 

K2 

whence 

/ A s = -fOOEK . (4.27) 

0k 


Operating on both sides with 

Dt 


DA' = -f(k)K 2 DE _ EK^ Df(k) + f(k)EK 2 Dk _ f(k)E D_ ( K 2) . 

Dt 0k Dt 0k Dt j} k 2 Dt 0k Dt (4.28) 


Noting that f(k) = Dk, and using (3.27a)-(3.27c) , 
Dt 3k Dt 


DA = -f(k)K 2 DE + [ -EK 2 3f + f(k)EK 2 ] Dk 
Dt 0k Dt 0k 3k g k 2 Dt 


-f (k)E 2 [kDk + A £* + I m £m .] 

0k Dt Dt ^ Dt 


= -f (k)K 2 DE + EK 2 f f - 3£ 1 . T -k 3u _ f 3v 1 

S Dt Ik - k 8k J 1 8? 8? 1 


+ f CK) E .2 • [k (k liL + l iJL) + Z (k j2i + £ Jbl ) 

0c 3x 9x 3y 9y 

+ JL m fk in. + 

^ 3z 



(4.29) 


Using (4.24a), (4.29) becomes 



7,-c\ [ 1a 2 * + m\ - [k (k is. + Jt 3*1 

Dt Bk J Bk 2 Bk 3x 3y 


+ £ fk iiL + £ ill + I m fk iH. + £ ill] 

3y 3y ^ 9z 9z 


EK 

Bk 


If 


“] 

3k 


f-k in. - t iy ] . 

3x 3x 


(4.30) 


Equation (4.30) can be simplified considerably by using the expressions for E 


given by (2.10) and for^by (4.27): 

& - 4 s Vo •c > + M[I-il][-k *»-£.£]. 


Dt 


0k k 3k 


3x 


3x 


Two choices of f(k) are of particular interest. 

1) f (k) = 1 . 

Then 'X' = ^ 1 ) = JL_ is wave action (Andrews and McIntyre, 1978), 


to 


^1) = - a 2 k 4 
4 Bk 


For this choice of f, (4.31) becomes 


(4.31) 


(4.32) 



(4.33) 


This important result (also given by Young and Rhines, 1980), states that even 
when the mean flow (u, v) depends on y and z, wave action is conserved locally , 
with the relevant flow of wave action occurring with the group velocity. 


2) f(k) = - Bk, whence 

'A' = ^(2) = EK 2 = P , 

so that 2) i s the wave potential enstrophy density given in (2.19). 
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From (4.31), 


+ v 3 • (tf 2 > • c) = o 

3t j 

- (4>34) 


This result, (again given by Young and Rhines, 1980) states that wave enstrophy 
is conserved locally for arbitrary (slow) variations of the mean flow (u, v ). 


- 28 - 




5 . Connection with Wave Activity, Potential Vorticity and the Eliassen-Palm 
Flux; Alternative Derivations 
5a. Eliassen-Palm Flux Vector 

-> 

The Eliassen-Palm flux vector F = (Fy, F z ) can be expressed (see the 
Appendix) as: 


F„ = e 


= p -z/h ~^r~W 
3y 3x 


F = e 


= -z/H 1 TrjT 


^ 3x 3z 


(5.1a) 


(5.1b) 


Here the primes denote perturbations from the zonal mean, and the overbars a 
zonal average. From our basic solution (2.10) for a Rossby wave, 


F y - lA 2 M , 

F = i -U 2 km , 

Z 'S' 2 

-> 

or F = 1a^ (kJl, _L km) . 

2 "s' 

Now define ^3) = p = _1_ K 4 = -k 'aO) , 

e 4fi 

and write down the intrinsic group velocity from (4 .12a )-( 4 . 12c) : 

^ _ ( r * 2f5k£ 2f5km\ 

K 4 1h< 4 


(5.2) 


(5.3) 


(5.4) 


From (5.2)-(5.4), 



or 


(C"A 


(3)) 


1 " 


-> 

F 


(5.5) 


Where 1 denotes a projection onto the y-z plane. Thus, the intrinsic group 
velocity times the wave enstrophy, divided by g, gives the Eliassen-Palm flux 

vector when projected into the y-z plane. Since only the y and z components 

-> -> 

enter (5.5), C > is identical to C if the basic state meridional velocity v vanishes. 

5b. Conservation of wave activity in the case of homogeneity in the x direction. 

When v = 0 and the wave amplitude and wavenumber (k, l, m) are independent 

of x , it is possible to derive a conservation law for wave activity directly from 

equation (2.9), without explicit consideration of the variation of A and (k, £, m) 

in y and z. We start by multiplying (2.9) by _L q * and zonally averging: 

B 

— ( — (q') 2 ) + v"q" = 0 . (5.6) 

at 28 

From (2.6b), 

v'q' = HI 8 2 iT + HI + 1 e z / H HI 1_ (e _z / R HI) 

3x 0 X 2 3 x 3y2 ^ 9x 3z 3z 

= H (HTh!) + I e z /H Hfe-z/n HThI) 

3y 3x 3y ^ 3z 3z 3x 


= (using (5.1a), (5.1b)) e z ^ H — F. + e z ^ H — F , 

3y y 3z z 

where we have used the fact thst ^ is assumed constent in z. 
Using (5.5), (5.7) becomes 


7^ = • (e z/H cV 3 )) , 


(5.7) 


(5.8) 


- 30 - 




-> -> _ -> 
where C can be used instead of because v = 0, and denotes the operator 

f j , JL.). Equation (5.8) relates the flux of potential vorticity to the divergence 
8y 3z 

of the group velocity times a quantity e 2 ^ A^) related to wave enstrophy. 

Now, from (2.11), our basic solution (2.9) gives 

_L_ (q') 2 = J_ K 4 A 2 e z / H = e z / H , 

2 6 4 6 

so we define *A^ 4 ) by 

^(4) = e z/H2(3) = e z/.H I P = i_ K 4 A 2 e z /H = 1_ (O? . 

8 46 26 (5.9) 


From (5.8) and (5.9), (5.6) becomes 


i-'A (4) + V, • (C tf*>) = 0 
3t 1 


(5.10) 


a( 4 ) is referred to by Edmon et al . (1980) as wave activity . Note that the 

local conservation law (5.10) only involves the two dimensional divergence 

V, = rather than the three dimensional form in, for example, (4.34). 

3y 3z -> 

This is because and C are both independent of x by assumption. It is 

interesting that this derivation does not involve any explicit calculation of 
-> ^ 

just how C and depend an y and z, and in particular how A and (k, £, m) 

vary in space. It is similar in this regard to a derivation of (4.34) presented 
by Hoskins et al. (1983). 
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APPENDIX 


1. Derivation of the quasigeostrophic potential vorticity equation in z 
coordinates. 

Our starting point is the quasigeostrophic equations in pressure coordinates, 
with constant static stablitity, from Lorenz (1960): 


JL£ = - J( 0) - a) 

3t 3p 


— (V^) = - j u>, - j (<|>, f) 

8t 

2 _> _> 
g V^z* = V 2 • (f V 2 <p) . 


-> -> 
V 2 • (f v D ) 


(A.l) 


(A. 2) 
(A. 3) 


In these equations, 0 is potential temperature, 0 S the reference state 

potential temperature (depending on p only), is the geostrophic streamfunction , 

“> 2 2 
V 2 the horizontal gradient operator and V 2 the horizontal Laplacian, V 2 ^ 

the vorticity, 00 = All, z* is altitude and f = 2 ft sin<f>, where Q is the angular 

dt 

rotation rate of the earth and 4> the latitude. J (A,B) represents the Jacobian 
of the variables A and B. Also, 


-> 

Vo 


-> 

V D 


+ = 0 
9P 


(A. 4) 


-> 

where vp is the divergent part of the wind. 

Neglecting the effects of spherical geometry, and neglecting the latitudinal 
variation of f except in the term J(\|>, f) in (A. 2), (A.1)-(A.3) become 

36 + 3ft 80 - 3ifr 39 = - a) 30s (A. 5) 

3t 3x 3y 3y 3x 3p 


0 ( v2\(;) + Ai ^2 ^ - At B Ai - f A^ = 0 

3t ^ 3 X gy 3y 3 X 3 X 0 3p (A. 6) 


gz* = f 0 ty . (A. 7) 
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In these equations y is the northward distance and x the eastward distance on a 
plane targent to the earth (the so-called "beta-plane") at latitude cf> 0 . (For 
a rigorous derivation of this approximation, see Chapter 6 of Pedlosky, 1979); 
f has been approximated in (A. 6) as f = f 0 + By, where f Q = sin<j> 0 . In 
(A. 7) f is taken simply as f Q . 

Introducing the new vertical coordinate (as in Section 2). 

z = - H to (p/p Q ), (A. 8) 

where H is the scale height, and p Q = 1000 mb, we have the identities 

JL=l z JL = :JiA_ = -H e z/H _a_ (k.9) 

9p 9p 9z p 9z p Q 9z 

and 

a)= dp = Jpdz = _3£ w=: -£ 0 e -z/H w (A. 10) 

dt 8z dt 3z H 

where w = is the vertical velocity in z coordinates, 
dt 


Using the equation of state pa = RT, where a is the volume per unit mass 

* 

and R the gas constant, and the hydrostatic balance = H2, where g is the 

3P g 

gravitational acceleration, we have 




Hence (A. 7) becomes 


T = Z£ f 

R ° 8p 


(A. 11) 


or 


T = Hfo 

3i|) 

R 

9z 


(A. 12) 


We convert (A. 5) to a prognostic equation for temperature by multiplying 

it by (— 2-) K , where k = R/ C , C being the specific heat at constant pressure: 
Po P P 
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(A. 13) 


3T + 3i|i 3T - 3ft 3T = - ^ (_p)< 99 fi 

3t 3x 3y 3y 3x p Q 3p 

Further, 

- a) (^.) K l§s = - 03 (JP) K: = - id jfe [(_£.)* <L!s ] . 

Po 9 P Po dp 9 P Po dz 


Defining (as in Section 2) 


S = (_£)< !®s, 

Po dz 


and using (A. 10), we have 


(A. 14) 


- 0 ) 


(_E) K lis = 

Po 9 P 


- w S 


so that (A. 13) becomes 


j5T+HjjT-j^j>T = _wS 

3t 3x 3y 3y 3x 


(A. 15) 


Finally, the term -f Q — in (A. 6) must be rewritten using (A.9) and (A. 10): 

3p 

_f 3u = _f 3 [3£ . W 1 . _ f 3z 3_ r3£ . „] . _ f e z/H J_ w j > 

3p 3p 3z 3p 3z 3z 3z (A. 16) 


whence (A. 6) becomes: 



(A. 17) 


The conservation law for (quasigeostrophic) potential vorticity is derived 
from (A. 12), (A. 15) and (A. 17): 


w = 1 3T _ 3\J> 3T + 3t|> 3T } 

5 3t 3x 3y 3y 3x 

which, using (A. 12) becomes 
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w 


(A. 18) 


= JLol {- a 2 4 - - 14 a 2 4 + .14 IjL] 

RS 3t 3z 3x 3z 3y 3y 3z 3x 


Corabining (A. 17) and (A. 18) gives, (again treating S as a constant): 


1_ ( v§ ip) + 14 1_ ( v?4) - 14 1- ( + b 14 

3t 3x 3y 3y 3x 3x 


- Lo^A e z / H 1_ { e -z / H [- d 2 4 - 14 1 2 _4_ + 14 IjL] } = o . (a.19) 

RS 3z 3z 3t 3x 3z 3y 3y 3z 3x 


We write the last term in (A.19) as a sum of three terms: 


f n ^H ( e z/H d_ | e ~z/H _ 3 ^_ 4 _} + e z / H 1 _ { e -z /H 14 4 } 

RS 3z 3z 3t 3z 3x 3z 3y 


- e z/H — { e _Zi/H 14 ill. }] = H [T, + T, + To] • (A. 20) 

3z 3y 3z3x RS 1 1 J 


T, - !- { e z/H 1- (e~ z/H 14) } , (A. 21a) 

1 3t 3z 3z 


= 24 e z / H JL (e“ z/H 92 4 } + 3 2 4 . 3 2 4 

3x 3z 3z3y 3z3y 3z3x 


= 14 !_ [e z / H 1_ fe _z / H 14) ] + g2 4 • IjL (A. 2 lb) 

3x 3y 3z 3z 3z 3y 3z 3x 


14 e z/H 1- {e _z/H LlL } - 114- . I 2 !. 

3y 3z 3z3x 3z3x 3z 3y 


= - 14 !_ [ e z /H 1_ (e-z/H 24)] - • jjfc (A. 21c) 

3y 3x 3z 3z 3z 3y 3z 3x 


Introducing the parameter S from equation (2,5) of the text, (A.19) through 
(A. 21c) may be combined: 
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JL [ V ? ip + I e z/H i_ (e" z / H ii)] 

3t z 'g' 3z 3z 

# 

+ M JL [v?ifr + I e z/H JL (e“ z/H 1^)] 

3y 3x ^ 3z 3z 

- li A- + I e z/H JL [ e -z/H M)]v+ 0 A£ = 0 . 

3y 3x ^ 3z 3z 3x 

This is just 


fi_ - it i. + it AJ q- « o 

3t 3y 3x 3x 3y 


with 


q = f + 3y + Vpi|> + i e z / H — (e Z ^ H A£) 

° z ^ 3z 3z 


We have thus derived equations (2.2) and (2.3) in the text. 


(A. 22) 


(A. 23) 


(A. 24) 
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2. Eliassen Palm Flux Vector in £n(p) coordinates 

The Eliassen-Palm flux vector, in pressure coordinates, is given by 
Edmon et al . (1980) for the beta plane: 



(A. 25) 


where primes denote wave quantities, the overbar a zonal average, and u, v are 

the geostrophic wind components: u = v = JNj, while 0 = T(P.o) K , where 

3y 3x p 

j n 

T is given by (A. 12). The static stability is given by a = — s , with 0 

d P s 

the reference state potential temperature. 

-> 

To obtain F in z coordinates (see equation (2.1)), we work with the (two- 

A> 

dimensional) divergence of F : 

Ky + Kp = + 9 _ (WTW) + ±_ (f vT) . 

3y 3p 3y 3x 3y 3p ° a 

Using the definition of S given in the text, and equation (A. 9), 

0 = dQ s = d6 q 3z = _ d9 s H e z/H = - s (_Po) K H_ e 2 /** 

dp dz 3p dz p 0 P Po (A. 26) 

whence 


f- (f 0 SE) - (-S- e z/K ) f- (£ 0 (£_)' • (- Po) e- z '») 

9p a Po 9z S Po H 

= e z / H JL (to e _z / H ) = (using (A. 12)) e z / H JL (i e“ z / H ) . 

9z S 3z ^ 3x 3z (A. 27) 

Using (A. 27), 

ity + itp = e z / H [JL (e" z / H jjtlll ) + JL (e _z / H i W' fjQ ] 

3y 3p 3y 3x 3y 3z 'g" 3x 3z 


= e z/H [t_ F + — F 1 . (A. 28) 

3y y 3z 2 

Equation (A. 28) defines (Fy, F z ) in the z coordinate system. They are identical 
to (5.1a) and (5.1b) of the text. 
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